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PHYLOGENY AND DISJUNCTION IN ROSCOEA 
(ZINGIBERACEAE) 


C. NGAMRIABSAKUL*, M. F. NEWMAN? & Q. C. B. CRONK*f 


A phylogenetic study of Roscoea (Zingiberaceae), a high-altitude genus of an otherwise 
tropical plant family, was undertaken using sequence data from the internal 
transcribed spacers (ITS) of the nuclear ribosomal DNA (nrDNA). Two species of 
Cautleya and two species of Curcuma were used as outgroups. This resulted in an 
aligned matrix of 436bp (ITSI, 203bp; ITS2, 233bp). Sequence divergence of ITS] and 
ITS2 within the ingroup ranged from 0-]3.9% and 0-7.6% respectively. The results 
suggest that Roscoea is monophyletic (BS=99%; DI= >3) with the genus Cautleya as 
sister group. Roscoea itself is divided into two sister clades which correlate with 
geography: a ‘Chinese’ clade (BS =67%; DI= +2) and a ‘Himalayan’ clade (BS = 59%; 
DI= +l). These two groups are disjunct across the ‘Brahmaputra gap’, a region in 
which no Roscoea spp. have been recorded. The only species which occurs on both 
sides of the Brahmaputra gap is Roscoea tibetica. However, the western populations of 
Roscoea tibetica (from Bhutan) show numerous morphological differences. It is 
therefore possible that Bhutanese R. tibetica represents a distinct taxon, possibly more 
closely allied to Himalayan species. 


Keywords. Biogeography, China, cladistics, ITS, the Himalaya. 


INTRODUCTION 


Roscoea is one of a group of five genera in Hedychieae (Zingiberaceae) which possess 
versatile anthers. The members of the group are Camptandra, Cautleya, Curcuma, 
Paracautleya, and Roscoea. They all occur in tropical regions or low-altitude sites, 
except the truly alpine genus, Roscoea. Cautleya and Roscoea have sometimes been 
confused by inexperienced observers. Indeed, these two genera occur in similar habi- 
tats, and have a similar habit with orchid-like flowers. Nevertheless, there are many 
characters separating these genera as pointed out by Cowley (I982), e.g. lateral 
petals are free from the claw of labellum in Roscoea while they are joined to the 
labellum for about half their length in Cautleya. Roscoea, the high-altitude genus of 
Zingiberaceae, comprises ]8 species (Cowley, I982; Cowley & Baker, I996). It occurs 
along the Himalaya from the west (Kashmir), to the east (south-western China), 
between |200 and 4880m (Cowley, ॥982). Roscoea grows in drier and cooler environ- 
ments than other Zingiberaceae, in places that are more exposed to extremes of 
climate. Unlike some other members of Zingiberaceae, Roscoea has closed leaf- 
sheaths (Spearing, I977). 

The internal transcribed spacer (ITS) sequences of nuclear ribosomal DNA are 
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well established as being useful in systematics (Baldwin, ॥992). ITS regions have 
rates of substitution that are useful for evaluating generic and species level relation- 
ships in plants (Baldwin et al., I995). Many investigations have been carried out 
using these regions, for instance in Asteraceae (Baldwin, I992), Apiaceae (Downie 
& Katz-Downie, I996), Gesneriaceae (Moller & Cronk, I997a, 997b), and Aralia- 
ceae (Wen et al., I998). In Zingiberaceae, there are currently phylogenetic studies 
going on at Royal Botanic Garden Edinburgh (RBGE) using ITS regions as a 
source of phylogenetic information. These regions have proved to be useful for 
studying the evolutionary relationships of the family at the species level (e.g. Curcuma 
spp., Ardiyani, I997; Alpinia spp., Rangsiruji, I999). 

This study aimed to confirm the monophyly of Roscoea, and the relationship 
between Roscoea and Cautleya, using living collections in the Royal Botanic Gardens 
(Edinburgh and Kew). It was hoped that, by combining data from ITS regions with 
distribution records and information on geological history, the study would give 
insights into the evolution of Roscoea and its sister genera: how a tropical plant 
family has colonized temperate regions. 


MATERIALS AND METHODS 
Ingroup taxa 


Eight species of Roscoea cultivated in the RBGE were verified by using the identifi- 
cation key and species descriptions of Cowley (॥982). Fresh leaf material of one 
plant representing each accession was taken for a total DNA extraction. Voucher 
specimens were prepared, flowers were also preserved in Kew cocktail (water 5.5 
units; methanol 3.5 units; glycerol 0.5 units) and both were deposited at the Royal 
Botanic Garden Edinburgh herbarium (E). DNA extracts of another eight species 
were taken from living plants at the Royal Botanic Gardens Kew (RBGK), and 
DNA aliquots were kindly provided by Dr Mark Chase (Table ]). These species 
represent all major areas of distribution of the genus (Fig. |). The remaining species 
that are no longer in cultivation are Roscoea nepalensis, Roscoea forrestii and Roscoea 
debilis. Attempts were made several times to acquire DNA of Roscoea nepalensis 
and Roscoea forrestii from dried herbarium specimens but unfortunately these failed. 
Roscoea nepalensis is an endemic species of central Nepal, near Jumla. It is thought 
that this species might be allied to those others from central Nepal, i.e. Roscoea 
capitata and Roscoea ganeshensis. On the other hand, Roscoea forrestii is one of the 
species that only occurs in south-central China. Although the DNA sample of 
Roscoea debilis from RBGK was thought to be genuine, it turned out to be a variant 
of Roscoea tibetica after closer examination. These ]5 species comprise most (83%) 
of the genus (total 8 species). 


Outgroup taxa 


There are currently 5 accepted names in Cautleya, though there may be fewer than 
five species (Kumar, ]994; Larsen et al., I998). C. carthcartii is probably just a 


TABLE ]. Accessions of Curcuma, Cautleya and Roscoea examined for ITS] and ITS2 sequence variation. * Number as shown in Fig. | the 
distribution map of Roscoea. PRBGE is Royal Botanic Garden Edinburgh; RBGK is Royal Botanic Gardens Kew. The distribution is given 
first and the locality of the plant sampled in this study is then given in brackets 


Distribution Royal Botanic Garden? Genbank 
accession 
(Number)* Taxon (locality sampled ) accession number number 
ITSI ITS2 

(I) Curcuma amada Roxb. SE Asia, India ( Kerala) RBGE ॥98] 000] AFI922I8 AFI922]9 
(2) Curcuma parviflora Wall. Burma, Thailand (Sukhothai) RBGE ॥985 ]66] AFI92220 AFI9222] 
(3) Cautleya gracilis (Sm.) Dandy India, Nepal, Bhutan, China, RBGE ॥982 0532 AFI92222 AFI92223 

Burma, Thailand (not known) 
(4) Cautleya spicata (Sm.) Baker India, Nepal, Bhutan, China, E0006739 AFI92224 AFI92225 

Burma (Nepal) (RBGE herbarium 

specimen) 

(5) Roscoea alpina Royle Nepal, Tibet, Bhutan, India RBGE 986 ]08 AFI92226 #७॥77१|92227 

(Himachal Pradesh) 
(6) Roscoea auriculata K. Schum. India, Nepal, Tibet (not known) RBGE ]969 9652 ॥]॥"॥92228 AFI92229 
(7) Roscoea australis Cowley Burma (Mount Victoria) RBGE ॥983 09]3 AFI92230 AFI9223] 
(8) Roscoea brandisii (Baker) K. Schum. India (Meghalaya) RBGK ]99/ 5649 AFI92232 AFI92233 
(9) Roscoea capitata Sm. Nepal (Ganesh Himal) RBGK ]992 2299 AFI92234 AFI92235 
(I0) Roscoea cautleoides Gagnep. China (Yunnan) RBGE 99] 0649 AFI92236 AFI92237 
(llI) Roscoea ganeshensis Cowley & W. J. Baker Nepal (Ganesh Himal) RBGK ॥992 2303 AFI92238 AFI92239 
(I2) Roscoea humeana Balf. f. & W. W. Sm. China (Yunnan) RBGE ॥985 60 AFI92240 AFI9224] 
(I3) Roscoea praecox K. Schum. China ( Yunnan) RBGK ॥994 35]] AFI92242 AFI92243 
(I4) Roscoea purpurea Royle India, Nepal, Bhutan (Ganesh RBGK ]992 230 AFI92244 AFI92245 

Himal) 
(]5) Roscoea schneideriana ( Loes.) Cowley China (Yunnan) RBGK ॥990 3345 ॥॥"॥92246 AFI92247 
(l6) Roscoea scillifolia (Gagnep.) Cowley China (not known) RBGE ॥979 4045 ॥॥"॥92248 AFI92249 
(l7) Roscoea tibetica Batalin Tibet, Bhutan, China (Yunnan) RBGE ॥985 I59 AFI92250 AFI9225] 
(I8) Roscoea tumjensis Cowley Nepal (Ganesh Himal) RBGK ॥992 230] AFI92252 AFI92253 
(I9) Rosocea wardiii Cowley India, Tibet, Burma (Yunnan) RBGE ॥987/ ]608 AFI92254 AFI92255 
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FIG. . Simplified geographical distribution of Roscoea species described to date (number 
5 to ]9 referring to the species listed in Table ]; number 20 = Roscoea debilis, number 2 = 
Roscoea forrestii, number 22 = Roscoea nepalensis). The position of the number is an indication 
of the species (note. two species, i.e. Roscoea alpina and Roscoea purpurea are widespread 
along the Himalaya). Cautleya gracilis (number 3) and Cautleya spicata (number 4) occur 
both in the Himalaya and China. Curcuma species numbers (] and 2) only indicate the origin 
of samples. Arrows show the course of the Brahmaputra river. 


robust form of C. gracilis while C. robusta is possibly synonymous with C. spicata 
(Smith, 994). The number of Curcuma spp. is less certain, partly because many 
species of Curcuma have long been widely cultivated, causing doubts on the justifi- 
cation of these species. Nonetheless, it is estimated at 50 species worldwide (Larsen 
et al., l998). Two species of Cautleya (C. gracilis and C. spicata) and two species of 
Curcuma (C. amada and C. parviflora) were chosen as outgroups because living 
collections of these plants are available at RBGE. As mentioned in the introduction, 
Cautleya is morphologically very similar to Roscoea. Its strong affinity with Roscoea 
necessitates further outgroup species which are distantly enough related to allow 
unequivocal rooting of the phylogenetic tree. Curcuma spp. were then chosen on the 
grounds that they possess versatile anthers, a shared distinct character of five genera 
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in Hedychieae, including Roscoea and Cautleya, but are clearly different from 
Roscoea and Cautleya in other characters. An attempt was made to obtain 
Paracautleya’s DNA, a monotypic genus from South India, from a dried herbarium 
specimen. Unfortunately, this was not successful. Fresh leaf material of Camptandra 
(4 species, Larsen et al., I998) from Malaysia (Ibrahim, pers. comm.) was not 
available, so it was not included In this study. 


DNA extraction 


Fresh leaf materials were kept in silica gel-filled plastic bags and stored at 0°C 
overnight in a refrigerator before extraction, to destarch the leaf tissue (starch may 
interfere with subsequent operations performed using the DNA). Total DNA extrac- 
tion was carried out using the modified CTAB procedure of Doyle & Doyle (I987) 
sometimes with further purification using the QIAquick® PCR purification kit 
(Qiagen Ltd, Dorking, Surrey, UK). All the samples of the study were obtained 
from fresh leaves, except Cautleya spicata which was taken from a dried herbarium 
specimen. 

Small scale total genomic DNA extraction using CTAB (Doyle & Doyle, 987) 
as a detergent gives lower levels of enzyme inhibition than other methods (Scott & 
Bendich, ]994). The modified protocol is as follows. 

A portion of leaf about lcm’ was cut into many small pieces, and put into a l.Sml 
microcentrifuge tube and about 50mg of purified sand and 200ul of 2x CTAB extrac- 
tion buffer were added. The leaf tissue was ground with a plastic pestle until a 
homogeneous slurry was formed. A further 800ul 2xCTAB was then added. The 
contents were mixed gently, and the tube was incubated at 65°C for 30min. The tube 
was allowed to cool to ambient temperature before adding 200ul of wet chloroform 
(chloroform 24 units; octan-l-ol | unit). 

The solution was mixed gently 4 or 5 times and centrifuged for 2min at ]3,000rpm. 
The aqueous (upper) phase was removed to a clean tube and re-extracted with 200ul 
wet chloroform. Again this was mixed gently to obtain a momentary single phase 
and centrifuged for 2min at ]3,000rpm. In another clean tube with the aqueous 
phase, 600ul cold (—20°C) propan-2-ol was added and the contents were mixed 
gently to precipitate the nucleic acids. After I0—lSmin at room temperature, the 
pellet of nucleic acids was precipitated by centrifuging for 2min at ]3,000rpm. The 
supernatant was removed and Iml of wash buffer (76% ethanol, IOmM ammonium 
acetate) was added. The tube was left for at least 30min to remove the 2x CTAB 
from the pellet. The supernatant was then aspirated as much as possible after the 
tube was centrifuged for 2min at I3,000rpm. Next, the pellet was dried completely 
by using an incubator drying oven for lOmin at 50°C. Lastly the pellet was dissolved 
in 30—SOul of sterile distilled water to obtain DNA concentration between I0—30ng/ul 
and stored at —20°C until required. 
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PCR amplification and sequencing strategy 


Double-stranded DNAs of the complete ITS regions In each genomic DNA were 
amplified by the polymerase chain reaction method (PCR ) using 2 primers, ITS 5P 
and ITS 8P (Moller & Cronk, ॥99/8). The primer sequences were (5’ to 3’), ITS 
S5P=GGA AGG AGA AGT CGT AAC AAG G and ITS 8P=CAC GCT TCT 
CCA GAC TAC A. The reaction (total volume=50ul) contained (in order of 
addition) 32.5ul of sterile distilled water, 5.0ul of l0x Dynazyme® reaction buffer 
(lx: IOmM Tris-HCl, pH 8.8 at 25°C, l.5SmM MegCl,, 50mM KCI], 0.I% Triton 
X-l00; Finnzymes Oy, Espoo, Finland), l.0ul of a mix of each dNTP at IOmM 
(final concentration 200M ) (Sigma Chemicals, Poole, Dorset, UK), 5.0ul of each 
primer at IOuM (final concentration luM) (Oswel DNA Service, Southampton, 
UK), a l.Oul aliquot of unquantified total genomic (template) DNA and 0.5ul (IU) 
of Dynazyme II thermostable DNA polymerase (Finnzymes Oy, Espoo, Finland). 
PCR amplification of the ITS region was carried out in 0.2ml microcentrifuge tubes 
in a Perkin Elmer thermal cycler. Each PCR reaction cycle proceeded as follows: ], 
Imin at 94°C to denature the double-stranded template DNA; 2, 2min at 55°C to 
anneal primers to single-stranded template DNA; and 3, Imin at 72°C to extend 
primers. The first cycle was preceded by an initial denaturation step of 3min at 94°C. 
Each set of reactions was monitored by the inclusion of a negative (no template 
DNA) control. Five microlitres of each double-stranded DNA PCR product were 
resolved by electrophoresis in ].5% agarose gel using Ix TBE as the gel buffer. 
Successful PCR resulted in a single band of ethidium bromide corporated-DNA 
viewed under ultraviolet (UV) light corresponding to approximately 700bp. The 
PCR product was then purified using the QIAquick® PCR purification kit. 

Purified PCR products were sequenced using a dye terminator cycle-sequencing 
ready-reaction kit (Perkin Elmer, Applied Biosystems Division, Warrington, UK ), 
with AmpliTag DNA polymerase, FS, according to the manufacturer's recommen- 
dation. Sequencing products were analysed on an ABI 377 Prism Automatic DNA 
Sequencer (Perkin Elmer, Applied Biosystems Division, Foster City, CA, USA), 
according to the manual supplied. Each reaction was 20ul in volume and contained 
(in order of addition) 6ul of sterile distilled water, Sul of Reaction Mix, lul of primer 
at 3.2uM and Sul of purified PCR product. For each taxon forward and reverse 
sequencing reactions were performed for sequence confirmation. Sequencing primers 
were ITS 5P, ITS 8P and in addition ITS 3P (Moller & Cronk, I997a) and a modifi- 
cation suitable for Zingiberaceae, ITS 2K (Rangsiruji, I999) were also used. All 
primers were synthesized by and purchased from Oswel DNA Service, Southampton, 
UK. The primer sequences were, ITS 3P=GCA TCG ATG AAG AAC GTA GC 
and ITS 2K =GGC ACA ACT TGC GTT CAA AG. 


Sequence analysis 


All sequences were verified by comparison of their forward and reverse sequences. 
Sequence boundaries of both internal transcribed spacers of all taxa were determined 
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by comparison with published rDNA sequence data for Daucus carota, Vicia faba 
(Yokota et al., I989) and Alpinia spp. (Rangsiruji, I999). ITS] and ITS2 of each 
species are deposited in GenBank (accession number AFI8695-AFI862I3, see 
Table ]). Both ITS regions were aligned using the CLUSTAL option in the multiple 
alignment program Sequence Navigator® Version ].0.] (Perkin Elmer, Applied 
Biosystems Division, Foster City, CA, USA), with minor manual adjustments. The 
G + C content was determined by inspection, and transition/transversion ratio using 
MacClade Version 3.0.] (Maddison & Maddison, I992). Sequence divergences 
among taxa were calculated using the DISTANCE MATRIX option in PAUP 
Version 3..] (Swofford, I993), based on unambiguously alignable regions (Fig. 2). 


Phylogenetic analysis 


Phylogenetic trees were generated using PAUP Version 3.].] (Swofford, I993), run 
on a Power Macintosh 6400/200 computer with character states unordered. The 
branch-and-bound search option, which guarantees to find the shortest tree or trees, 
was selected with MulTrees and furthest addition sequence options. 

Bootstrap analyses (Felsenstein, I985) were performed using PAUP, set to branch- 
and-bound search option and ]000 replicates. Decay indices (DI: Bremer, ]988; 
Donoghue et al., I992) for individual clades were obtained by comparing the strict 
consensus of all equal-length trees up to four steps longer than the shortest tree, 
using branch-and-bound search option. Descriptive statistics reflecting the amount 
of phylogenetic signal in the parsimony analyses were given by the consistency index 
(CI: Kluge & Farris, ]969), retention index (RI: Farris, I989), and the resulting 
rescaled consistency index (RC: Swofford, ॥993). Additionally, the g, statistics (Hillis 
& Huelsenbeck, ]I992) were obtained by calculating the tree-length distribution of 
]0,000 random trees using Random Trees under PAUP to assess the amount of 
phylogenetic signal in the data set, in comparison to random noise. 

For all analyses of sequence data, gaps (indels) were treated as missing data (Soltis 
& Kuzoff, I995; Susanna et al., I995; Downie & Katz-Downie, I996). Indels were 
scored as a separate presence/absence character and added to the sequence data 
matrix (Wojciechowski ef al., I993; Oxelman & Lidén, ॥995). To investigate the 
effect of these additional data, a separate analysis without indels scored as characters 
was undertaken. Character state changes were weighted equally, except for one analy- 
sis in which character-state weighted parsimony was implemented: transversions were 
weighted over transitions by a factor of ].7, corresponding to an average of the 
transition/transversion ratio of ITS] and ITS2. 


Mapping the distribution area of Roscoea 


Two hundred and eighty eight records of locations of Roscoea spp. were taken from 
all the herbarium sheets at E, including extra locations taken from a revision of 
Roscoea (Cowley, I982). These data were entered into PANDORA (a taxonomic 
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Curcuma amada TTTGGGATGAGTCCTCCAGAGACCCTGTGTGATGATTGCGGAGTCGCGTGAAAGCGCCGCGTCAATCA---TTTGC [430] 
Curcuma parviflora TTTGGGATGAGTCCTCAAGAGACCCTATGTGAT---TGCAGAGTCGGACGAAAGCGATGTGTCAATCATCATTTGC [49] 
Cautleya gracilis TTTGGGAAT-GTCCTCAAGAGACCCTGTGTGAT-~ ~-TGTGATGTCGTGTGAAAGTGCCGTGTCCATCA--AATTIGT [45] 
Cautleya spicata TTTGGGAAT-GTCCTCAAGAGACCCTGTGTGAT— --TGTGATGTCGTGTGAAAGTGCCGTGTCCATCA- -AATTGT [45] 
R.cautleoides TITAGGATT-GTCCTCAAGAGACCCCGTGTGAT - - -TGTGACGTCGTGCGAAAGTGCCGTGTCCATCA- ~AATTGT [4759] 
R.wardii TTTAGGATT-GTCCTCAAGAGACCCCGTGTGAT---—TGTGACGTCGTGCGAAAGTGCCGTGTCCATCA--AATTGT [425] 
R.humeana TITAGGATT-GTCCTCAAGAGACCCCGTGTGAT-— —--TGTGACGTCGTGCGAAAGTGCCGTGTCCATCA- -AATTGT [45] 
R.praecox TTTAGGATT-GTCCTCAAGAGACCCCGTGTGAT- —--TGTGACGTCGTGCGAAAGTGCCGTGTCCATCA- -AATTGT [45] 
R.australis TTTAGGATT-GTCCTCAAGAGACCCCGTGTGAT - - -CGTGACGTCGTGCGAAAGTGCCGCGTCCATCA- -AATTGT [45] 
R.scillifolia TTTAGGATT-GTCCTCAAGAGACCCCGTGTGAT-— —- -TGCGACGTCGTGCGAAAGTGCCGCGTCCATCA- -AATTIGT [43] 
R.schneideriana TTTAGGATT-GTCCTCAAGAGACCCCGTGTGAT- - -TGTGATGTCGTGCGAAAGTGCCGTGTCCATCA- -AATTGT [425] 
R.tibetica TTTAGGATT-GTCCTCAAGAGACCCCGTGTGAT-— - -CGTGATGTGGTGCGAAAGTGCCGTGTCCATCA--ATTTGT [4i3] 
R.capitata TTTAGGATT--TCCTCAAGAGACCCTGTGTGAT---TGTGATATCGTGCGAAAGTGCCGTGTCCATCA--AATTGT [4]2 |] 
R.tumjensis TTTAGGAT'I—-'TCC'ICAAGAGACCCCC'TG'IGAT--- TGTGATATCGTGTGAAAGTGCCGTGTCCATCA--AATTGT [42} 
R.ganeshensis ATTACGATT-—-TCCTCAAGAGACCCTGTGTGAT- —--TGTGATATCGTGTGAAAGTGCCGTGTCCATCA~-AATTGT [43] 
R.purpurea TTTACGATT-—-TCCTCAAGAGACCCCGTGTGAT-— --TGTGATGTGGTGTGAAAGTGCCGTGTCCATCA- -AATTGT {43 |] 
R.brandisii TTTACGATT— -TCCTCAAGAGACCCCGTGTGAT- —--TGTGATGTCGTGTGAAAGTGCCGTGTCCATCA--AATTGT {42] 
R.auriculata TTTAGGATT-——TCCTCAAGAGACCCCGTGTGAT-— —- -TGTGATGCGGTGTGAAAGCCCCGTGTCCATCA--AATTGT [4]2] 
R.alpina TITAGGATT--TCCTCAAGAGACCCCGTGTGAT---TGTGATGTCGTGCGAAAGTGCCGTGTCCATCA--AATTGT [4]2] 


FIG. 2E. Sequence data matrix of aligned ITS] and ITS2 regions of nuclear ribosomal DNA of ]9 taxa of Zingiberaceae. Nucleotide sequences 
are displayed from 5’ to 3’. ITS] ranges from site | to 203 and ITS2 ranges from site 204 to 436. Uncertain nucleotide states are coded according 
to PAUP conventions (Swofford, I993): n=A/C/T/G, k=G/T, r=A/G, s=C/G, w=A/T, y=C/T, m=A/C; hyphens denote alignment gaps; 
numbers in italic print above the nucleotide matrix, ranging from | to ]4, indicate the number and position of alignment gaps; numbers in 
square brackets at the end of sequences indicate the actual spacer length of the combined region of ITS] plus ITS2. 
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database system by Richard Pankhurst and Martin Pullan, RBGE: 
www.rbge.org.uk/research/pandora.home) at RBGE. The latitude-longitude format 
data in PANDORA were then exported and modified for use with MapPad (a 
freeware program by John Keltner, NOAA Paleoclimatology Program: 
http://www.ngdc.noaa.gov/paleo/softlib.html. The simplified distribution map of 
Roscoea is shown का Fig. | (including the outgroups). 


RESULTS 
Sequence analysis 


Alignment of internal transcribed spacer sequences of the ]9 taxa analysed resulted 
in a 436bp long data matrix (Fig. 2); its characteristics (including G+C content) 
are given in Table 2. The number of unresolved bases ranged from 0 to 5bp 
per sequence. 

The lengths of ITS] and ITS2 were, on average, I89.7 and 225.lbp. Alignment of 
all taxa required the insertion of ]4 gaps of | to 5bp length, 8 in ITS] and 6 in ITS2 


TABLE 2. Sequence characteristics of ITS] and ITS2 regions of ]9 taxa of Zingiberaceae 


Parameter ITS] ITS2 ITSI and ITS2 
Length range (total) (bp) l88—200 224-230 4] 2--430 
Length mean (total) (bp) l89.74 225.05 4[4.79 
Length range (ingroup) (bp) ]88-]90 224-225 4]2-4]5 
Length mean (ingroup) (bp) ]88.93 224.53 43.47 
Length range (outgroup) (bp) 90-200 225-230 4]5--430 
Length mean (outgroup) (bp) ]92.75 22] 4]9.75 
Aligned length (bp) 203 233 436 

G +C content range (%) 47.34--55.79 53.07—59.56 5l.55—57.35 
G +C content mean (%) 52.43 56.64 54.73 
Sequence divergence (ingroup) (%) 0—3.86 0-7.58 0-9./5 
Sequence divergence (in/outgroup) (%) 3.2l-9.22 4.462] 4.58-]8.47 
Number of indels (ingroup) 3 l 4 

Number of indels (total) 8 6 l4 

Size of indel (ingroup) |-2 I l-2 

Size of indel (total) l-5 l-3 l-5 
Number of variable sites (%) 67(33) 73(3.33) ॥40(32.0) 
Number of constant sites (%) ]36(67) ]60 (68.67) 296( 67.90) 
Number of informative site (%) 27(3.30) 43(8.45) 70(6.05) 
Number of autapomorphic sites (%) 40(9.70) 30(2.88) 70(6.05) 
Transitions (minimum) 50 7] ]2] 
Tranversions (minimum) 40 32 de 
Transitions/tranversions L.25 2.2] l.68 
Skewness of tree-length distribution 

(g, value for 0,000 random trees) --.022 -].063 --].509 
Average number of steps per character 0.448 0.446 0.447 
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of which 5 and 2, respectively, were potentially informative. The lengths of aligned 
ITS! and ITS2 regions were 203 and 233bp respectively. Of these aligned sites, 296 
(67.90%) were constant, 70 (l6.05%) had at least two nucleotide states In two or 
more sequences and were potentially informative phylogenetically, and 70 (I6.05%) 
were autapomorphies (Table 2). 

Sequence divergence of ITS] and ITS2 between ingroups ranged from 0-3.9% 
and from 0-7.6% respectively. Sequence divergence between ingroups and outgroups 
showed that ITS2 was marginally more variable at 4.5—2].0% than ITS] at 
3.2-]9.2%. Pairwise comparison of individual taxa across both spacer regions 
revealed 0—9.7% sequence divergence within the ingroup and 4.6--]8.4% divergence 
between ingroup and outgroup taxa analysed (Table 2). The maximum sequence 
variation between Roscoea accessions was 9.7% (40 character changes) between R. 
praecox and R. ganeshensis. Sequences of R. cautleoides, R. wardii and R. humeana 
were identical. 


Phylogenetic analysis 


Parsimony analysis of aligned ITS sequences using equally weighted character states 
yielded five most parsimonious trees when coded indels were added to the data 
matrix. The strict consensus tree was computed (Fig. 3), with 2]3 steps when all 
uninformative characters were included, ॥36 steps with autapomorphies excluded, 
with CIs of 0.8]2 and 0.706, respectively. These were higher than the expected empiri- 
cal value of 0.559 calculated from 60 phylogenetic studies for ]9 taxa (Sanderson & 
Donoghue, ॥989). The RI was 0.793, and thus the RC was 0.644 with, and 0.560 
without, uninformative characters. 

The average number of nucleotide substitutions per character was low, with 0.447 
indicating a low saturation of base substitution. The homoplasy index (HI) of the 
present data matrix was low (HI=0.88). 

Thirty-three character changes separated the Cautleya/Roscoea clade from 
Curcuma spp. (BS=00, DI= >3). The ingroup Roscoea spp. was separated from 
Cautleya spp. by nine character changes (one indel) (BS=99, DI= >3). Roscoea 
spp. formed 2 distinct groups: l, a Chinese clade comprising eight species from China 
and Burma, separated by two character changes (BS=67, DI=2); 2, a Himalayan 
clade with seven species from the Himalaya, separated by four character changes 
(one indel ) (BS=59, DI=]) (Fig. 3). In the Chinese clade, the relationship of species 
was fairly well resolved, with bootstrap values ranging from 59 to 75% and decay 
index values of | to 2. However, the relationship of a terminal branch in this Chinese 
group was unresolved due to the lack of sequence variation, forming a four species 
polytomy (R. cautleoides, R. wardii, R. humeana, R. praecox) separated from R. 
australis by two character changes (BS =75, DI=]). The Himalayan clade contained 
two subclades: l, R. capitata, R. tumjensis and R. ganeshensis with seven character 
changes (BS=70, DI=]); 2, R. auriculata and R. alpina by one character change 
(BS =53, DI=]) with R. purpurea and R. brandisii unresolved. Of the seven poten- 
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FIG. 3. Strict consensus tree based on five most parsimonious trees for ]5 Roscoea, two 
Cautleya and two Curcuma taxa of 2]3 steps length based on parsimony analysis of the 
combined ITS] and ITS2 sequence data plus the coded indels. Upper numbers are bootstrap 
values of 000 replicates. Lower (boldface) numbers are decay indices (CI =0.8l2; RI=0.793; 


RC=0.644). 


tially informative indels, four were congruent with the tree topology of the strict 


CONSENSUS tree. 


Exclusion of the coded indels from the combined ITS] and ITS2 data matrix 
resulted in eight most parsimonious trees of ]95 steps (l25 steps excluding unin- 
formative characters; CI=0.8]5; RI=0.787; RC=0.642). The strict consensus tree 
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differed from the strict consensus tree obtained with the addition of coded indels 
only in the collapse of the Himalayan clade grouping all seven species from the 
Himalaya while the two subclades within remained. 

The transition/transversion ratio was ].25 for ITS] and 2.2 for ITS2, and ].68 
for the combined data matrix. Altering the character weights to ].7:] to accommo- 
date this ratio and reanalysing the data (coded indels excluded) in a parsimony 
analysis gave a single most parsimonious tree (Fig. 4). 


DISCUSSION 
Molecular evolution of ITS in Roscoea 


The internal transcribed spacers of ribosomal DNA of the Roscoea spp. investigated 
have evolved mainly by base substitution. Only four indels occurred in the DNA 
studied, of |—2 bases in length. The levels of sequence variation among the Roscoea 
spp. are similar to those infrageneric levels found in other angiosperms. Sequence 
divergence within the Roscoea spp. ranged from 0 to 3.8% for ITS] and 0 to 7.5% 
for ITS2. In species of Saintpaulia (Gesneriaceae) for example, the range of sequence 
divergence was from 0 to I7.6% for ITS] and 0 to 3.9% for ITS2 (Moller & Cronk, 
997b) and in species of Alpinia (Zingiberaceae), sequence variation ranged from 0 
to 20.9% for ITS] and 0 to I9.7% for ITS2 (Rangsiruji, 999). However, a group 
of four species remained unresolved because the level of sequence divergence was 
too low for unequivocal phylogenetic resolution. Indeed, three of these species had 
identical sequences. Other similar studies have such unresolved groups (Kim et al., 
996; Moller & Cronk, I997b) and this is generally attributed to rapid radiation, 
especially on islands or in newly created ecological niches. The Chinese and 
Himalayan mountains apparently represent a recent range extension for the predomi- 
nantly tropical family Zingiberaceae. This may have induced processes of adaptive 
radiation similar to those found on islands. These regions have been affected by the 
continuous uplift of the Himalaya since the collision of the Indian and Asian plates 
०.52 to 45.8Ma B.P. (Rowley, I998). 

Although the spacers show considerable variation at higher levels of the taxonomic 
hierarchy, they are thought to be important in post-transcriptional processing, and 
thus are conserved to some extent (Liu & Schardl, ॥994; Van Nues et al., I994). It 
is interesting to find that speciation (as in the Saintpaulia ionantha complex of Möller 
& Cronk, I997b and in the Roscoea cautleoides complex in this study) has been able 
to outstrip variation in the comparatively fast-evolving ITS region. Moller & Cronk 
(997b) have suggested that, where divergence times are short compared to rDNA 
homogenization rates, ITS variation will appear highly conservative; on the other 
hand, where divergence times are long compared to rDNA homogenization rates, 
ITS variation will appear disproportionately variable. 
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FIG. 4. The single most parsimonious tree obtained from the weighting of transitions and 
tranversions (one of the five trees found in the unweighted search). Numbers above branches 
indicate number of character changes shared amongst taxa (branch length, from unweighted 
analysis), including autapomorphic changes. Bars and numbers associated indicate the indels 
and their positions In the sequences (see Fig. 2). 


Roscoea and Cautleya 


The phylogenetic trees resulting from this study show that Roscoea is monophyletic 
and Cautleya is its sister group. This is supported by a preliminary phylogenetic 
study of Hedychieae which included R. cautleoides, R. purpurea, C. gracilis and 
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twenty six species from the other ten genera of Hedychieae (Searle & Hedderson, in 
press). The relatively low sequence divergence of species of Cautleya from those of 
Roscoea suggests a close relationship between the two genera. This is supported by 
their similar morphology and overlapping distribution area. Roscoea spp. can be 
found between 200 to 4880m and Cautleya spp. (Kumar, ॥994; Wu & Larsen, in 
press) between 900 to 3]00॥7 above sea level. However, Cowley (I982) pointed out 
that some clear distinguishing features exist. Roscoea has no true petiole, its lateral 
petals are free from the claw of the labellum and it has an elongate capsule, while 
Cautleya has a true petiole, the lateral petals are joined to the claw of the labellum 
and it has a round capsule. In addition, Roscoea has small (c.3 x lcm), fusiform, 
fascicled tuber roots, whereas in Cautleya, the tuber roots are cylindrical. The closed 
leaf sheath (Spearing, I977) of all Roscoea spp. and Cautleya gracilis also suggests 
a relationship. 

It is also interesting to consider which of the Chinese and Himalayan groups of 
Roscoea is more closely related to Cautleya. Morphologically, Cautleya spp. are 
superficially similar to Roscoea cautleoides. However, the phylogenetic analyses pre- 
sented here suggest that neither the Chinese nor the Himalayan clade can be con- 
sidered as more closely related to Cautleya, as they are sister groups. 


Two groups in Roscoea 


The strict consensus tree of five most parsimonious trees resulting from the combined 
ITS sequences and coded indels clearly shows not only that Roscoea is monophyletic, 
but that it is also divided into two distinct groups. Seven species from China and 
one species from Burma form the first group (Chinese clade) (Fig. 3), while the rest 
(seven species from the Himalaya) form the second group (Himalayan clade). These 
two groups are supported by morphological characters as shown in Table 3. In order 
to explain this divergence, we need to examine the distribution of Cautleya, the sister 
group of Roscoea. Cautleya is not only found with Roscoea at lower levels of the 
Himalaya and in south-central China, but is also recorded from high-altitude sites 
of nearby tropical mountains, in Burma and in the north of Thailand (Larsen, I980). 
However, the geographical centre of the present distribution of Roscoea and Cautleya 
is Assam (including all provinces in north-eastern India; level 3 botanical region, 
Hollis & Brummit, ]992), as shown in Fig. |. Assam is also the centre of diversity 
of the related genus Hedychium. There are 39 species of Hedychium in India, of 
which 35 occur in Assam (Jain & Prakash, I995). It is possible that Roscoea orig- 
inated in Assam, and spread east and west along the nearest mountain ranges, thus 
accounting for the separate Chinese and Himalayan groups. This is supported by a 
single maximally likely tree showing that a clade of Roscoea/Cautleya shares an 
ancestor with Hedychium species clade (Searle & Hedderson, in press). Smitinand 
et al. (I970) reported that Anaphalis margaritacea ( L.) Benth. & Hook. f. ssp. margar- 
itacea is mainly distributed in the cool temperate zone of eastern Asia (including the 
Himalaya), and in North America. On finding the species in northern Thailand, he 
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TABLE 3. The distinguishing characters between the two groups of Roscoea spp. 


Chinese group Himalayan group 

l. Sheathing leaf number 3-5 Sheathing leaf number 0-2 

2. Leaf number 0-4 Leaf number >4 

3. Leaf base not auriculate, except R. Leaf base auriculate, except R. capitata 

australis, R. tibetica 

4. Leaves almost forming rosette, except Leaves not forming rosette 

R. cautleoides 

5. Corolla tube length <6cm, except R. Corolla tube length >6cm 

humeana, R. australis 

6. Appendage tip obtuse (R. australis, R. Appendage tip pointed 

scillifolia: obtuse-pointed, R. schneideriana: 

globular) 

7. Epigynous gland length <5mm Epigynous gland length >5mm, except R. 
auriculata, R. capitata 

8. Seed aril deeply lacerate, except R. Seed aril shallowly lacerate 

scillifolia, R. australis, R. wardii 

9. Ratio labellum length/dorsal petal Ratio labellum length/dorsal petal length 

length <], except R. schneideriana, R. >] 


wardii, R. cautleoides 


suggested that the species may have spread southwards along the high mountains 
of the Indo-China Peninsula to Thailand and Vietnam. Similar migration along 
mountain dispersal routes may have occurred in Roscoea and Cautleya. 

All of the species in the Chinese clade (except R. australis, Burma) are found in 
Yunnan province (mostly in Lijing and Dali) and some extend to parts of Sichuan. 
The data suggest that this is an area of rapid evolution of a complex of Roscoea 
spp. On the other hand, the area of greatest diversity of the Himalayan clade is in 
central Nepal. One particular area is Ganesh Himal (Cowley & Baker, I996) which 
accounts for up to five species among eight species in the entire Himalayan region. 
These data give an indication of the priority of land protection and preservation for 
the authorities concerned. 


The Brahmaputra gap 


The distribution of Roscoea (Fig. ]) is strikingly discontinuous. There are no records 
from that part of Assam where the Brahmaputra river flows south around the eastern 
end of the Himalayan chain. Interestingly, this gap in the distribution coincides with 
the boundary between the Chinese and Himalayan clades. Although it is possible 
that the Brahmaputra gap is an artefact of undercollection, it is also possible that 
it represents a genuine phytogeographical boundary. 

The region of the Brahmaputra gap is known to be undercollected, as the area 
has been historically inaccessible. Rao (994) suggests that 30% of north-eastern 
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India (not including Arunachal Pradesh) has been only casually surveyed. More 
collecting in this region is therefore badly needed. 

It is also possible that this area really has no Roscoea spp. Although the Himalayan 
mountains form a continuous, geologically connected chain, here the eastern 
Himalaya rise rather abruptly from the plain without a distinct sub-Himalayan zone 
(Rao, I994). This abrupt rise of the mountain range and its horseshoe shape may 
serve as a barrier between the two sides of the area. Thus the disjunct distribution 
of Roscoea, between two sides of north-eastern India, may be genuine along with 
other examples of Indian disjunctions (Rao, I994), Le Nymphaea pygmaea Ait. (with 
Siberia, N China), J/licium cambodiana Hance (with Southern Indo-China), 
Mitrastemon yamamotoi Makino (with Japan, Sumatra) and Dendrobium bensoniae 
Reichb. (with Burma, Thailand). 


Roscoea tibetica as a transgressor species 


As mentioned above, the species of Roscoea fall into two groups, with either an 
eastern or a western distribution. The only exception to this is Roscoea tibetica which 
occurs in both the eastern (China and south-eastern Tibet) and the western area 
(Bhutan and in nearby Tibet) (Fig. 5). There are two possible explanations for this: 
l, that R. tibetica is a genuine transgressor which crosses the phytogeographical 


S Brahmaputra river 





FIG. 5. Distribution map of Roscoea tibetica showing the discontinuity between Chinese 
and Bhutanese populations. Arrows show the course of the Brahmaputra river. 
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TABLE 4. The distinguishing characters between the two geographically distinct populations 


of Roscoea tibetica 


Chinese populations 


l. Leaf number 22 

2. Calyx longer than bract 

3. Corolla tube long, exserted from calyx 
4. Labellum shorter than lateral petals 

5. Labellum usually divided more than half 
6. Labellum drying dark purple or pink (in 
herbarium specimens) 

7. Labellum throat with white lines 

8. Lateral petal tip acute 

9. Appendage tip obtuse 

0. Stigma with long hairs at tip 


Bhutanese populations 


Leaf number >3 

Bract and calyx equal 

Corolla tube short, within calyx 
Labellum longer than lateral petals 
Labellum divided less than half 
Labellum drying purple (in herbarium 
specimens) 

Labellum throat without white lines 
Lateral petal tip obtuse 

Appendage tip pointed 

Stigma with short hairs at tip 


boundary; or 2, the Bhutanese populations of R. tibetica represent a separate species, 
possibly more closely allied to Himalayan ones. 

The accession used in this study comes from China and groups with the Chinese 
clade. It would be very interesting to obtain material from Bhutan. There is some 
evidence of morphological difference between the Bhutanese and Chinese specimens 
(Table 4). Further studies of Bhutanese Roscoea tibetica remain a priority. 
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